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ABSTRACT 



o 
o 

I We study correlations amongst tidal fields originated by the large-scale distribution of 

matter in the Universe. The two-point tidal correlation is described as a rank-4 tensor, 
\ whose elements can be written in terms of four fundamental scalar functions ranging, 

with respect the spatial separation, from purely transversal to purely longitudinal cor- 
relations. Tidal fields, both on galaxy and cluster scales, reveal to be correlated over 
distances larger than the mass-density correlation lenght, though traceless tidal fields 
show anti-correlation between diagonal terms along orthogonal directions. The cross- 
■ correlation between mass and tidal field is also analyzed. These results are relevant 

CNl ' for galaxy formation and the interpretation of large scale weak lensing phenomena. 

OO ' 

o 
o 

O ■ 1 INTRODUCTION 



Key words: galaxies: statistics - large-scale structure of Universe. 



The role played by cosmic tidal fields in the evolution of mass density fluctuations in the Universe has been recognized since 
the very first studies of gravitational instability dynamics (see e.g. Peebles 1980 for a historical perspective). Torques induced 



O ' by tidal fields presumably prompted the process of angular momentum acquisition by a proto-galactic halo (Hoyle 1949; 

White 1984; Heavens & Peacock 1988; Catelan & Theuns 1996). Tidal fields influenced the dynamics of the collapse around 
^ ' primordial dark matter fluctuations (Doroshkevich 1970; van de Weygaert & Babul 1994; van de Weygaert & Bertschinger 
1996; Bond & Myers 1996). Consequently, they also determined the mass and spatial distributions of virialized halos (Monaco 
1995; Catelan et al. 1998; Lee & Shandarin 1998; Jing 1998, 1999; Porciani, Catelan & Lacey 1999; Sheth & Tormen 1999). 
Only recently, however, the importance of spatial correlations amongst cosmic tidal fields has been evidenced. For instance, 
5_l ' our ability of reconstructing maps of the cosmic mass distribution from weak gravitational lensing measurements (Bacon, 
Refregier & Ellis 2000; Kaiser, Wilson & Luppino 2000; van Waerbeke et al. 2000; Wittman et al. 2000) is hampered by 
our ignorance about the distribution and alignment of galaxy shapes. As a first guess, one can assume that shapes and 
orientations of background galaxies are uncorrelated, and that the detected correlations amongst galaxy ellipticities are solely 
due to the foreground mass distribution. However, in the last few months, a series of observations and numerical simulations 
evidenced that both galaxy shapes (Croft & Metzler 2000; Brown et al. 2000) and angular momenta are intrinsically correlated 
(Heavens, Refregier & Heymans 2000; Pen, Lee & Seljak 2000). Theoretical models for shape alignments, whether induced by 
spin correlations or not, ascribe them to spatial correlations of the underlying tidal field (Catelan, Kamionkowski & Blandford 
2000; Crittenden et al. 2000). 

In this Letter we analyze two-point statistics of the cosmic tidal fields, addressing some fundamental issues that have 
been only marginally discussed in the aforementioned literature. We determine, for example, the size of a tidally correlated 
region, and the characteristic length of the cross-correlation between linear density fluctuations and tidal fields. If intrinsic 
alignments of galaxy shapes are ultimately determined by the cosmic tidal fields, these are the fundamental scales involved in 
the process. In §2, we briefly remind the definition of the cosmic tidal field. In §3 we compute the rank-4 tidal-tidal correlation 
tensor. Section 4 discusses astrophysical implications and contains our conclusions. 



2 COSMIC TIDAL FIELDS 

Let us consider the gravitational potential </?(<!) originated at the comoving Lagrangian position q by the primordial mass 
density fluctuations S in the Universe. We define ip through the Poisson equation, \7 2 ip = S, where the density S can be thought 
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as the primordial density field extrapolated to the present time. For our purposes, it is irrelevant whether tp is a Gaussian 
field or not. The tidal field T> a p is defined as a rank-2 symmetric tensor 

V a0 = dadpip = -J ^jp a Pfi^(p)e' p ' q , (1) 

where d a = d/dq a , p is the comoving Lagrangian wavevector and tp is the Fourier transform of tp. More rigorously, 2? a/ 3 
is known as deformation tensor; the tidal field is properly the traceless rank-2 tensor defined in terms of the gravitational 
potential tp as T a p = [dadp — (l/3)5 a V 2 ]</?. In a similar fashion, it is possible to introduce a third tensor describing the 
velocity shear. In fact, the flow of an irrotational fluid as the cosmic dust can be described by introducing a peculiar velocity 
potential, ip v . The shear tensor is then defined as [d a d/3 — (l/3)5 a/ 3V 2 ]v3„. During the linear regime we consider in this paper, 
(p v — tp, and we are left with only one tensor d a dptp which can be called, indifferently, tidal or shear tensor, if we are not 
interested exclusively in its traceless part. 



3 CORRELATIONS AMONGST TIDAL FIELDS 

Let us now compute explicitly the tidal-tidal correlation, which turns out to be a rank-4 symmetric tensor: 

£ a /3 7CT (q) = (f a /3(qi)T' 7CT (q2)> = d a dp <9 7 <9 CT (, v {q) , (2) 

where q = qi— q2|, and £, v is the gravitational potential correlation function, £ v (q) = (¥>(qi)<^(q2)) = (2-7r)~ 3 J dp P v (p) e lpq . 
Here, Ptp(p) is the gravitational potential power spectrum, (tpCp)tp(p')) = (27r) 3 fo(p + p')P<p(p), and Sd denotes the Dirac 
function. The explicit calculation of the right hand side of eq.(|2j) is somewhat toilsome; the reader should bear in mind that 
d a = (dq/dq a )(d/dq) = (q a /q)d/dq, then apply iteratively this identity to £ v and its radial derivatives. At the end of the day, 
it turns out that the tidal-tidal correlation tensor can be decomposed as follows 

£ a /3 7CT (q) = T(q) T a p 1<7 + V(q) V a f3 1C r + 5(g) S a p 1<7 , (3) 
where, respectively, 

T( 9 )=C'(9)-?C(9) + S^(9)-5^(9). n?) = ^C(3)-l^(9) + 4C(«). «(«) = k(«)-k(«); ( 4 ) 
q q A q A q q z q A q A q a 

with f L = d£,p (q) /dq and so forth. In addition, 

T q a qi3 5 7 5<t 

Va/3~<o = 5— 7(T "I ^— 0^ CT H — 7/ 3 -I 5-ia/3 H 5— 0-, a H 5— QCT (Ol 

q z q 2, q z q z q z q z 

SaPja = SapSja + 5 la 5f3cr + 5a<r6jf3 , (7) 

where 8 a p is the Kronecker symbol. ^ Crittenden et al. (2000) report [their eq.(34)] an alternative expression for eq.(^). 
The equivalence of the two results can be shown using text-book relations like, for example, f^(<?) — — — 2q ,_1 ^(g), 
and i' v {q) = -q- X i v {q) + (2n 2 )- 1 J dpp^Ptfp) cos(qp)/qp, and = {q/5)J 5 (q), and C'(<?) = £«(«) - (4/5) J 8 (5) and 

so forth. Here &(q) = (<S(qi)<5(q 2 )), and J n (q) = ng~ n J* dco lu"- 1 &(u;). Eventually, one gets: T = & + (5/2) J 3 - (7/2) J 5 , 
V = (1/2) J 8 - (1/2)J 3 , S = (1/6)J 3 - (1/10) J B . 

The most striking feature appearing in eq.(^) is the anisotropy of the tidal-tidal correlation function: the spatial derivatives 
combine in different ways depending on the orientation of the separation vector with respect to the frame of reference. The 
anisotropic structure of £ a ^ 70 - correlations can be evidenced by explicitly writing its 81 components. For simplicity, we choose 
the reference frame in such a way that the positive z-axis lies along the separation vector q = qi — q 2 . In this way, the only 
non- vanishing components are 

£1111 (?) = £2222(5) =3 5( g ) = \jz(q) - ^J 5 (q) , (8) 

£1122(9) = £2211(5) = £1212(5) = S(g) = ij 3 (g) - -jJ B (g) , (9) 

£1133(5) = £2233(5) = £1313(9) = £2323(5) = S(q) + V(q) = ~J 3 {q) - |j S (<z) , (10) 

£3333(5) = 35(5) + 6V(q) + T(q) = &(<?) - I J B (q) , (11) 

5 

* Note that, in linear perturbation theory, the tidal correlation tensor evolves proportionally to D(t) 2 , with D{t) the growth factor of 
density perturbations. 
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(the components of the tensor are invariant under permutations of the indices). In total, there are 21 non-zero elements 
in the whole tensor. With respect to the separation-axis, we could refer eq.(^) to eq.(|ll|) as the 'orthogonal 1-orthogonal 
1', 'orthogonal 1-orthogonal 2', 'parallel-orthogonal' and 'parallel-parallel' tidal correlations: they can be considered as the 
basic building blocks of the more general tensor in eq.(^|). These fundamental functions are plotted in the left panel of Fig. [j] 
considering a ACDM cosmology. Q 

Recently, the problem of estimating the size of the tidally correlated sphere about any given point has become of relevance. 
To this purpose, one has somehow to consistently get rid of the anisotropies of the tidal correlation tensor. The simplest way 
to attain this is to average the correlation tensor over all the directions, {C a/ 3 7tT )47r = [(1/15)T + (2/3) V + S] 5 a| g 7C r = 
(1/15) £s Safs-ya, and the contributions from the off-diagonal elements cancel out. Alternatively, one can saturate the free 
indices in £a,/3 70 -(q) in order to construct scalar (i.e. rotationally invariant) quantities. Obviously, there is not a unique way 
to do this. A first attempt might be to consider the simplest £ aaaa = ^aC^i unfortunately, this operation also eliminates 
purely tidal contributions from cross-terms like 9 2 9f ^ an d so forth. The next simplest scalar one might look at is <L aa pp\ but, 
trivially, this is the mass-mass correlation function, deprived, again, of any information about the actual tidal correlations 
(induced by the off diagonal elements of the tidal tensor). One is forced to consider higher-order scalars. The simplest of them 
is obtained saturating C a( 3 7(T with itself: £ Q( 3 7CT £ a /3 7( T ■ [Note that this quantity will be quadratic in the underlying correlations, 
so one might prefer instead to consider the scalar (£ q1 /3 7<t £q/3 7 o-) 1 ' /2 -] Notably, this is the same quantity to which Lee & Pen 
(2000) and Pen, Lee and Seljak (2000) ascribe the leading-order spin-spin correlations. 

Summarizing, we assume here that the scalar (C a 7 o-l£a/3 7 <T/81) 1/ ' 2 is an estimate of the typical tidal correlation about 
any given point in the Universe. We find, after saturating the indices, 

£ a/ 3 7CT (q)e: a/ 3 7CT (q) = (© a , a (qi)© 7CT (q 2 )> 2 = &(qf - | Jis +2 Jz{qf -4 J 3 (q) Js(q) + y Mqf ■ (12) 

A discrepancy with the results reported by Lee & Pen (2000) appears here. Their eq.(H5), in fact, claims that £ 2 l( 3 7(T = £ 2 , 
where the several other terms appearing in the right hand side of our eq.(jl^), all of the same order of magnitude, are neglected. 
(Rather, the correct expression would be C 2 ^^ = £f .) 

Another interesting quantity to study is the cross-correlation between mass concentrations and tidal fields. This is, in 
fact, a special kind of tidal-tidal correlation, since the mass fluctuations correspond to the trace of the tidal matrix. Mass- 
tidal correlations could be detected in the vicinity of galaxy clusters, e.g. through weak gravitational lensing. From the 
Poisson equation, the mass-tidal rank-2 correlation tensor can be obtained from the more general rank-4 tensor <t a p Jt7 via the 
saturation of two indices, for example 7 and a; we obtain 

£ a/377 (q) = (O a/3 (q 1 )C 77 (q 2 )> = [T(q) + 7V(q)\ ^ + [V(q) + 5S(q)} S a0 = Mq) - J 3 (q)] ^ + ± J 3 (q) 6 a0 . (13) 



3 



Again, a rotationally invariant scalar can be constructed from C Q( 3 77 (q) by squaring it; 



£c/3 77 (q)£ Q/ 3<T CT (q) = £<5(<7) 2 - - J3(q) £s(q) + - -h{.qf '■ (14) 

This function is plotted in the right panel of Fig.ji]. Note that [£a/3 77 (q)£ a /3 CT (T(q)/9] 1 '' 2 is an estimate of the typical mass-tidal 
cross-correlation. It is also of interest to consider the correlation of the traceless tidal field T a /3 7(T (q) = (J r a/3(qi)J r 7 <T(q2)): 

2 a( g 7(T (q) = £ a /3 7 ,T(q) — -5 7 ct ^^^^(q) — — 5a/3 £i7»7 7 o-(q) + -(S a/ 3 5 7 ct £i/i-?7?7(q) 



^4q) + ^[^(g)-^(g)](^5 7CT + ^^)+i[6(?)-2J 3 (g)]^5 7CT . (15) 



For q = (0,0, q), only the following components do not vanish, 

Xiiiifo) = T 2222 (g) = + ^J 3 (q) - Aj 5 ( g ) , (16) 

%U22(q) = W( 5 ) = -&(<?) - -^Mq) - ±Mq) , (17) 

T 1133 (?) = T 2233 (g) = - p 3 (q) + jU(q) , (18) 

^3333(9) = ^s(q) + p 3 (q) - p 5 (q) , (19) 
and also Xi 2 i 2 (g) = £i 2 i 2 (g), Ti 3 i 3 (g) = £1313 (?), < J 2 3 2 3(?) = £ 2 3 2 3(q). The corresponding correlation scalar is then 



t The model is specified by the following parameters. Mass density parameter O = 0.3, cosmological constant density parameter Q\ = 0.7, 
Hubble constant Ho = 100 fekm/s/Mpc, with h = 0.7. A cold dark matter (CDM) power spectrum of density fluctuations with primordial 
spectral index n = 1 is assumed. This power spectrum is normalized by requiring that rms density fluctuation in a 8 h~ 1 Mpc sphere is 
= 0.9 when linearly extrapolated to present-day. 
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Figure 1. Linear two-point correlations of the deformation tensor vs. spatial separation, in a ACDM cosmology. The separation vector 
lies along the z-axis. All the amplitudes are linearly extrapolated to the present time. Left panel. The solid line shows the correlation 
between components along the two directions perpendicular to the separation vector (01-02), £1122 = £1212 = £2211; the dotted line 
shows the correlation between diagonal components orthogonal to the separation vector (Ol-Ol), £1111 = £2222- The short-dashed line 
shows the correlation £1133 = £1313 = £3311 = £2233 etc. (P-O). The correlation of the diagonal component along the separation is 
plotted with a long-dashed line (P-P). Right panel. The correlation scalar (£ a 7 CT£a/3 7 o-) 1/ ' 2 (solid line) is compared to \£s\ (dotted 
line), and to the mass-tidal correlation (£ a( 3 77 £ a( 3 (T<T ) 1 ^ 2 (short dashed line). The typical tidal-tidal correlation, (£ a( 3 7tT £ a ^ 7(T /81) 1 / 2 , 
(long-dashed line) and mass-tidal correlation, (£ Q /3 77 £ a /3 (T < T /9) 1 / 2 , (dot-dashed line) are also plotted. 



4 2 

T Q/ 3 7CT (q)T a/ 3 7CT (q) = -6(g) +8 



7^3(9) - 7 J s(<?) 



fc(ff) + ^Mlf ~ 4J 3 (g) J B (q) + ^JrMY 
9 5 



(20) 



This quantity is plotted in the right panel of Fig. 



4 DISCUSSION AND CONCLUSIONS 

A tidal correlation length for (proto) dark matter haloes can be defined by smoothing the fluctuation field on the relevant 
scales, and determining the separations at which the normalized correlation functions, e.g. £a/3 7CT (q)/|£a/3 7 (T(0)| (the covariance 
functions divided by the corresponding variances), assume a chosen value C. At zero separation, we obtain £ a /3 7<T (0) = 
<t 2 /15 S a /3 J(T , where a 2 is the variance of density fluctuations. The correlation lengths q 1 / 2 and qi/ 3 (corresponding to C = 1/2 
and C — 1/3 respectively) are reported in Table 1 for a smoothing radius of 1 /i -1 Mpc (galaxy scales) and of 10 /i" 1 Mpc 
(cluster scales). Top-hat smoothing is adopted in every case. Note that, by definition, q x /2 and gi/ 3 keep constant during the 
linear evolution of the density field. For a ACDM cosmology, we found that, with respect to the separation, purely orthogonal 
tidal correlations (Ol-Ol and 01-02) tend to be stronger than the longitudinal tidal correlations (P-O and P-P); this holds 
true for the traceless tidal field as well. The latter shows systematic anti-correlations, being the 01-02 and P-O correlations 
negative on the relevant scales. In addition, the tidally correlated sphere around any given point, computed here via the scalar 
£ 2 ( 3 70 ., is bigger than the mass correlated sphere. 

Lee & Pen (2000) and Pen, Lee & Seljak (2000) used the expression £ 2 , ( 3 7<T = £ 2 to quantify the leading-order three- 
dimensional correlation amongst galaxy spin directions (|L ■ L'j 2 ) which, in turn, would model correlations amongst intrinsic 
galaxy shape alignments: a fundamental issue for reconstruction of dark matter maps from weak gravitational lensing mea- 
surements. A more comprehensive discussion of the implications of our eq.(|l2^) for the correlations amongst galaxy angular 
momenta and for the reconstruction of dark matter maps from weak gravitational lensing measurements is beyond the scope 
of this Letter. However, if intrinsic alignments of galaxy shapes are ultimately determined by the cosmic tidal fields, the scales 
reported in Table 1 must be the fundamental scales involved in the process. 

The results reported here do not assume that the gravitational potential is Gaussian distributed, but if it is so, then our 
analysis gives a complete description of the tidal-tidal statistics. 

CP acknowledges the support of the EC RTN network "The Physics of the Intergalactic Medium". PC thanks Marc 
Kamionkowski and Rupert Croft for discussions. This project has been partially funded by NSF AST-0096023. 
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Figure 2. As in Fig. |l| but for the traceless tidal field correlation T a ^ 7(T . Left panel. All symbols are like in Fig. (lj). The solid line 
represents T1122 = "-£2211; ^1212 has not been plotted since it coincides with solid-line curve in FigJlJ Analogous consk erations are valid 
for T1133 = X2233 (short-dashed line) and T1313 = X2323 (short-dashed line in Fig. [l]). The functions T1122 and T1133 are negative, 
indicating anti-correlation bewteen diagonal components of the tidal tensor. Right panel. The scalars ( t S Q ^ 7(T T C( ^ 70 -) 1 / 2 (solid line) and 
C£a/37CT^a/37<r/81) V 2 (long-dashed line) are compared to \£g\ (dotted line). 



Table 1. Correlation length (in h~ 1 Mpc) of the tidal field, for different smoothing 
scales, i?TH (same units). Negative numbers denote anti-correlation lengths. £ 2 = 
£ a /3 7(T £ a /37(T! M-T (mass-tidal)=£ a) 3 77 £ a) 3 (T(T . A ACDM cosmology is assumed. 
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